Short periodic orbit approach to resonances and the fractal Weyl law 
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We investigate the properties of the semiclassical short periodic orbit approach for the study of 
open quantum maps that was recently introduced in [M. Novaes, J.M. Pedrosa, D. Wisniacki, G.G. 
Carlo, and J.P. Keating, Phys. Rev. E 80, 035202(R) 2009]. We provide conclusive numerical 
evidence, for the paradigmatic systems of the open baker and cat maps, that by using this approach 
the dimensionality of the eigenvalue problem is reduced according to the fractal Weyl law. The 
method also reproduces the projectors |'0n)('0n |, which involves the right and left states associated 
with a given eigenvalue and is supported on the classical phase space repeller. 

PACS numbers: 05.45.Pq 



I. INTRODUCTION 

One of the cornerstones of the semiclassical approach 
to quantum mechanics is the Gutzwiller trace formula, 
which relates the fluctuating part of the quantum den- 
sity of states to the classical periodic orbits of a chaotic 
system [1]. However, this relation requires infinitely long 
orbits and involves divergent sums. More recently, pe- 
riodic orbits have been used to study quantum spectra 
in a different way, by constructing some special quantum 
states, called scar functions, which are adapted to the 
system's classical dynamics and provide a suitable basis 
for diagonalizing the Hamiltonian. This formulation has 
the advantage of using only a small number of short pe- 
riodic orbits, and shedding light on the phenomenon of 
scarring which is an anomalous localization of sta- 
tionary wave functions along periodic orbits. 

In quantum scattering a prominent role is played by 
resonances or quasibound states. These are eigenfunc- 
tions of the system with complex energy, whose (nega- 
tive) imaginary part is interpreted as a decay rate. In 
chaotic systems the number of states with a prescribed 
decay rate grows as a power of the energy which is con- 
jectured to be related to a fractal dimension of the clas- 
sical repeller, the set of initial conditions which remains 
trapped in the scattering region for all, positive and neg- 
ative, times. This fractal Weyl law has been investigated 
in several systems One the other hand, the (right) 

eigenfunctions are supported by the unstable manifold of 
this repeller p!ol, flU , and its scarring properties have also 
been under investigation 

EMI. 

The scattering analog of the short periodic orbit ap- 
proach has been introduced in [isj]. Scar functions were 
constructed which are concentrated on the periodic orbits 
but also extend along the unstable manifolds by means 
of a dynamical evolution up to the order of the system's 
Ehrenfest time. This basis is adequate for the calculation 
of the small fraction of resonances which have small decay 
rates, usually the most important ones, without having to 
consider the multitude of rapidly decaying states. Also, 
it provides a tool for studying localization on classical 



structures more directly, and allows a natural semiclassi- 
cal approach to resonance wave functions. 

This method was used in [15] to reproduce the main 
resonances of an open baker map for a specific value of 
h. In the present work we focus on the scaling of the 
method and show that it is compatible with the frac- 
tal Weyl law, not only for the baker but also for a cat 
map, which is more generic. We also discuss the ability 
of method to reproduce eigenfunctions. We find that a 
mixed quantity involving both types of functions (right 
and left), which has been introduced in [16], can be ac- 
curately reproduced. The paper is organized as follows. 
Section [TTl is focused on reviewing the quantum and the 
classical versions of the systems that we have used in our 
study: the open baker and cat maps. In Section [IIll we 
discuss the method and its results. Finally, we present 
our conclusions in Section HVl 



II. OPEN QUANTUM MAPS 

Maps are paradigmatic systems in classical and quan- 
tum chaos because of their simplicity |l7l4l9(. We will 
consider maps defined on the torus. When quantizing 
them, boundary conditions must be imposed for both the 
position and momentum representations. This amounts 
to taking (g + l|V^) = e''^''^^ {q\^l)) , and (p+l|V^) = 
e*^^^p(p|V^), with Xq-, Xp ^ [O^l)- This implies a Hilbert 
space of finite dimension N = (27rfi)~^, and the semi- 
classical limit is approached for large N. The system's 
propagator becomes a N x N matrix. The discrete 
set of position and momentum eigenstates is given by 
\Qj) = \U + Xq)/N) ^nd \pj) = \{j^Xp)/N) with 
j G {0,...,A^ — 1}. They are related by a discrete Fourier 
transform, i.e. {pk\qj) = -^e-^^^^^'+^^^^^+^^'^Z^ = 

We consider a certain region of the torus to be the 
opening, through which particles can escape. Its quan- 
tization is implemented by means of a projection op- 
erator P on its complement. We always choose a fi- 
nite strip parallel to the p axis and corresponding to a 



2 



range of q values, so that the projector P is quite sim- 
ple in position representation. If U is the propagator 
for the closed system, then U = PUP represents the 
open one. It will have N right eigenvectors |^^) and N 
left ones (^^|, which are orthogonal among themselves 
(^^1^^) = dnm- Their norm is arbitrary, but one may 
choose (^^1^^) = (^^1^^). 
The classical (tri)baker map 

(3(7, p/3) if < g < 1/3 

B{q^p) = { (3g - 1, {p + l)/3) if 1/3 < g < 2/3 (1) 
(3(7-2,(p + 2)/3) if 2/3 < g < 1 

is an area-preserving, uniformly hyperbolic, piecewise- 
linear and invertible map with Lyapunov exponent A = 
In 3. Following [20",'2l], the quantum version is defined in 
terms of the discrete Fourier transform in position repre- 
sentation as 
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(2) 



where antiperiodic boundary conditions are imposed, 
Xq = Xp = 1/2. For this system we always take the 
opening as the region 1/3 < q < 2/3. 
The classical cat maps are of the form 
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FIG. 1: Exact spectrum (circles) and our results (crosses), for 
the baker map. The pairs (N,Ns) consisting of the Hilbert 
space dimension N and the number of scar functions Ns are: 
a)(81, 51); b) (177, 105); c), (597, 231) and d) (1821, 471). We 
have chosen Ns ~ 4:N^^'^ where d is the dimension of the 
classical repeller. 
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where the Cij 



(3) 

must be integers to ensure continuity and 
the conditions TrC > 2 and detC = 1 are imposed to 
make the map hyperbolic and area-preserving. Here, we 
consider 
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3 2 



(4) 



for which the Lyapunov exponent is log (2 + a/3) and the 
stable and unstable directions are s = (—a/3, 1) and u = 
(\/3, 1). Quantization of cat maps was first introduced in 
Ref. [l3] and discussed in [22]. For the case considered 
here, with periodic boundary conditions = Xp = 0? 
this results in 




2i7r(Q2-Q'Q+Q'2)/Ar 



(5) 



where q = Q/N and q' = Q'/N. 

We note that for the baker map the opening corre- 
sponds to a cell in the Markov partition. As a con- 
sequence, the repeller is given in terms of an exactly 
self-similar fractal, the well known middle-third Cantor 
set. In that respect the cat map is more generic, since 
its stable and unstable manifolds intersect the opening 
transversally. 



III. METHOD AND RESULTS 

Scar functions are special wavefunctions constructed 
by taking into account classical information in the neigh- 
borhood of a periodic orbit [23'-^295. They have been de- 
veloped for closed systems and are the building blocks of 
the semiclassical theory of short periodic orbits, by means 
of which one can find eigenvalues and eigenfunctions of 
a quantum system starting from purely classical quanti- 
ties. For open systems they were introduced in Ref. [l5j. 
In this section we review their construction and show a 
few examples. 

Let 7 be a periodic orbit of an open map (it must 
therefore belong to the repeller) of fundamental period 
L, i.e. it consists of L different points in the torus: 

(<?o,Po), (<7i,Pi), ....((7l-i,Pl-i), (<7l,Pl) = (<?o,Po). (6) 

We associate with 7 a total of L scar functions. Initially, 
we define coherent states \qj,Pj) for each point of the 
orbit and a linear combination of them called a periodic 
orbit mode. 



1 

^ exp{-27ri(j4 - Nej)}\qj,pj). 



(7) 



i=o 



Here k G {0, . . . , i — 1} and 9j = J2i=o '^h where Si is 
the action acquired by the ^th coherent state in one step 
of the map. The total action of the orbit is = S'-y and 
= {NSj + k)/L. 
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FIG. 2: Exact spectrum (circles) and our results (crosses), 
for the cat map. Here, the pairs {N,Ns) are a)(114, 90); b) 
(415,198) and c), (1751,476). 



The right and left scar functions associated with the 
periodic orbit are defined through the propagation of 
these modes under the open map. Namely, 



7,/c 



and 



7,fel 



1 
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1^'), (8) 
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The constants A/"^'^ are chosen such that {i^^^k\^^,k) — 

(V^7,/elV^7,/e) (V^7,/c I V^7,/e) = ^' ^hc COSiuC is USCd tO 

introduce a smooth cutoff. The time scale of the propa- 
gation, r, is taken proportional to the system's Ehrenfest 
time. 

We use these functions to construct an approximate 
basis in the Hilbert space to diagonalize our propagator. 
We select a number of short periodic orbits that approxi- 
mately cover the repeller, and eigenvalues and eigenfunc- 
tions are obtained by solving a generalized eigenvalue 
problem [l5|. In this way, we isolate the relevant in- 
formation needed to construct just the long lived reso- 
nances, without calculating the others. According to the 
fractal Weyl law, the number of such resonances grows 
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FIG. 3: Husimi representation of exact eigenstates and the 
results of our method for the baker map ai N = 243 and 
\z\ = 0.895; right eigenstate in panels a) (exact) and b) (our 
method), while c) and (d) are analogous for the left eigenstate. 
The overlaps between these pairs are KV^iLl^fc)!^ = 0.727 and 
\{ipex\4^sc)\'^ — 0.556. Differences are underlined by means of 
magenta (dark gray) rectangles. 



(d) 







FIG. 4: Husimi representation of exact eigenstates and the 
results of our method for the cat map. Differences are un- 
derlined by green rectangles. Right eigenstate in panels a) 
(exact) and b) (our method), while c) and (d) are analogous 
for the left eigenstate. Here N = 100 and \z\ = 0.722. Over- 
laps are | (V^il |0fc> I' = 0.954 and | (V^i; |0?c> I' = 0.948. 



like 7V^/2 [3-9], where d is a fractal dimension of the 
classical repeller. Our method takes advantage of this 
fact: the number of scar functions that we need to ob- 
tain a reasonable approximation to the long-lived sector 
of the spectrum, denoted Ns, is of the order of N^^'^. 
The dimension of the matrix to be diagonalized is thus 
substantially reduced. 

Exact spectra and the results of our method are shown 
in Fig. [T]for the baker map. We notice a very clear gap 
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FIG. 5: Mixed representation hn{q,p) of the same eigenfunc- 
tions shown in Fig. [3]and|4l coming from: a) exact calcula- 
tion; b) our method. Results for the baker map are on the 
left and for the cat map are on the right. The agreement is 
much better than in Fig. OandHl for both systems. 

developing in the spectrum, which is reproduced by the 
method. This gap has been observed before [5], but the 
reasons for its existence have not yet been understood. 
The number of scar functions used was taken to scale 
with dimension as the fractal Weyl law: Ng ~ A^^/^, 
where d is the dimension of the classical repeller. The 
quality of the eigenvalues obtained does not deteriorate 
with indicating that this choice is correct and the 
long-lived sector of the spectrum indeed has a reduced 
effective dimensionality. 

Analogous results are presented for the cat map in Fig. 
[21 The quality of the individual eigenvalues is not as good 
as in the baker map. However, the scaling is preserved, 
i.e. the results confirm that we only need about N^^'^ scar 
functions for each value of N in order to reproduce the 
portion of the spectrum lying closest to the unit circle. 

It has been shown [ijj that the right eigenstates are lo- 
cated in the unstable manifold of the repeller, in the sense 
that their Husimi representations ^ are 
nearly zero outside that set. On the other hand, left 
eigenstates are supported by the stable manifold. Our 
method is based on functions that are approximately sup- 
ported by the repeller. Therefore, perhaps surprisingly. 



although it is able to provide very accurate eigenvalues, 
the corresponding eigenfunctions are not necessarily well 
reproduced. We show in Fig. [3] the right and left eigen- 
functions of the eigenvalue z = 0.895 of the baker map 
for N = 243. We see that the results from our method 
have significant differences compared to the exact ones, 
in regions away from the repeller. We also show analo- 
gous results for the cat map {z = 0.722, N = 100) in Fig. 
m Differences are not as noticeable in this case. 

The method is, however, able to accurately reproduce 
the mixed representation hn{q^p) = 
which was recently introduced in [l6|. The reason for 
that is that this phase space quantity is supported on the 
intersection of the individual supports, which is precisely 
the repeller where our scar functions live. In Fig. Owe 
can clearly see that the values of hn{q^p) computed with 
our method indeed coincide with the exact ones for both 
systems. 



IV. CONCLUSIONS 

We have shown in two paradigmatic models that our 
recently developed method [15] to obtain quantum reso- 
nances of chaotic systems from their classical properties 
indeed shows the same scaling as the fractal Weyl law. 
Namely, only a fraction N^^'^ (where d is the dimension 
of the classical repeller) of our special basis states are 
required to provide a good approximation to the long- 
lived sector of the spectrum. When the Hilbert space di- 
mension N becomes large, this actually represents only 
a small fraction of N. We found that the method gives 
good but not excellent approximations to the right and 
left resonant eigenstates independently, because of the 
very limited support of the basis states we use. However, 
a mixed quantity involving both right and left eigenstates 
can be successfully obtained. Perhaps with some modifi- 
cation the method could also reproduce them, but further 
investigation is needed. 
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